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Surface states in Tl/Si(lll)-(lx 1) and /3-Bi/Si(lll)-(\/3 x \/3) show non-Rashba-type spin split¬ 
ting. We study spin-transport properties in these surface states. First we construct tight-binding 
Hamiltonians for Tl/Si and Bi/Si surfaces, which respect crystallographic symmetries. As a result, 
we find specific terms in the Tl/Si surface Hamiltonian responsible for non-Rashba spin splitting. 
Using this model we calculate current-induced spin polarization in the Tl/Si Hamiltonian in order 
to see the effect of non-Rashba spin-orbit interaction. We found that the induced spin polarization 
is in-plane and perpendicular to the current, which is consequently the same with Rashba systems. 
We find that it follows from crystallographic symmetries. Furthermore, we numerically find bound 
states at the junction between two surface regions which have different signs of the spin-orbit inter¬ 
action parameters in the Bi/Si system and in the Tl/Si system. We explain these numerical results 
with the results of our analytical calculations. 
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I. INTRODUCTION 

The spin-orbit interaction (SOI) in crystals has various 
forms depending on materials and their symmetries. In 
a surface of a metal, one of the typical forms of the SOI 
is the Rashba SOI (a x [I| and is observed in metal 
surfaces such as in Au(lll) surface [2[. Here the Fermi 
surfaces consist of two concentric circles, and spins are 
tangential to the Fermi surfaces. This spin splittin^ffect 
has been observed on clean noble metal surfaces [3|-[6| and 
heavy group V elements M- On the other hand, bulk 
inversion asymmetry gives rise to another type of the SOI 
called Dresselhaus SOI [ 13 . The Dresselhaus SOI in two 
dimensions has the form: kx(Jx — kyCfy to the linear order 
in k. Its energy band splitting is similar to that of the 
Rashba SOI, but the spin direction is unlike the Rashba 
SOI, as has been observed in the GaAs(llO) surface 0- 

Since there are various types of SOI terms depend¬ 
ing on systems, we can expect rich physics from their 
interplay. For example, two-dimensional (2D) systems 
including both the Rashba and Dresselhaus SOIs with 
equal magnitude have an anomalously enhanced spin life¬ 
time 0- Because of nested Fermi surfaces with oppo¬ 
site spins, the spin life time is proposed to be largely 
enhanced at the nesting wave vector. This phenomenon 
was observed in 2D electron gas [TsI in semiconductor 
quantum wells la and is called persistent spin helix. 

Thanks to the low symmetries of material surfaces, 
they sometimes allow even other types of SOI. For exam¬ 
ple, non-Rashba-type surface states have been measured 
in Tl/Si(lll)-(1 xl) [l5| and /3-Bi/Si(lll)-(\/3x v^) sur¬ 
faces [13 by angle-resolved photoemission spectroscopy 
(ARPES). On the Tl/Si surface the spin-split states at 
the K point have the spins normal to the surface. In ad¬ 
dition, on the Bi/Si surface there is a “peculiar” Rashba 
splitting at the K points. It is called peculiar because 


the K points are not time-reversal invariant, whereas the 
conventional Rashba splitting appears only around time- 
reversal invariant k points 00- Furthermore, at the 
M point in the Bi/Si surface, the spin texture is not of 
Rashba type, but is similar to that of Dresselhaus type. 
Such unconventional systems may have new spin proper¬ 
ties which are yet to be discovered. 

In this paper, we theoretically explore new proper- 
ties due to the SOI in the Tl/Si(lll)-(lxl) and 
/3-Bi/Si(lll)-(v/3 X y/S) [T^ surfaces. First, we con¬ 
struct effective tight-binding Hamiltonians of the Tl/Si 
and Bi/Si surfaces. We verify qualitative agreement for 
energy bands and spin texture, between the experimental 
results and our results. Second, we explore spin proper¬ 
ties in the non-Rashba-type system, Tl/Si, such as the 
current-induced spin polarization and the persistent spin 
helix. We show that a charge current induces spin polar¬ 
ization in the system with non-Rashba-type spin split¬ 
ting. Although the Hamiltonian contains out-of-plane 
components, the induced spin is shown to be in-plane. 
In addition, when the Fermi energy is controlled in the 
Tl/Si surface so that there are two nested carrier pockets 
with opposite spins, we theoretically show that spin helix 
states can be realized at a nesting wave vector. 

Additionally, we study a junction between two surface 
regions which have different signs of the SOI parameters 
in the Bi/Si and the Tl/Si systems. It is motivated by 
the related work in topological insulators (TIs) [T^ . The 
junction of TIs with different sizes of the SOI is shown 
to exhibit a novel refraction phenomenon, and in addi¬ 
tion it was shown that there exist topologically-protected 
gapless interface states between two TIs when the SOI of 
the two TIs have opposite signs [T^. Because TIs and 
Rashba systems have the same form of the SOI, we ex¬ 
pect a similar behavior in Rashba systems. Motivated 
by this work, we study a junction between two surface 
regions with non-Rashba SOI of opposite signs. We see 
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that in some cases the junctions support bound states, 
with its spin directions different from the bulk states. 


II. CONSTRUCTION OF EFFECTIVE MODELS 
FOR THE NON-RASHBA SYSTEMS 

In this section we construct effective tight-binding 
models for two types of non-Rashba systems. These mod¬ 
els are intended to be minimal models sharing the same 
symmetry properties as the original systems. Hence they 
do not necessarily reproduce the band structures of the 
original materials quantitatively. 


A. Tl/Si(lll) surface 

We construct an effective tight-binding Hamiltonian 
for the Tl/Si surface [Fig.^a-l).]. We consider the model 
on the triangular lattice, which represents the T1 atoms 
as shown in Fig. [T]^a-2). We retain only the nearest- 
neighbor hoppings, and the x and y directions are defined 
as in Fig.[Ha-2). By taking into account the symmetries 
of the Tl/Si crystal, such as three-fold rotational symme¬ 
try Csz^ time-reversal symmetry, and mirror symmetry 



FIG. 1. (Color online) (a-1) Schematic illustration of the 
Si(lll)-(1 X1) surface with adsorption of one monolayer of 
Tl. The orange balls indicate the topmost Si atoms and the 
blue balls indicate the Tl atoms, (b-1) Schematic illustration 
of the Si(lll)-(A/3 X \/3) surface with adsorption of one mono- 
layer of Bi. The purple balls indicate the Bi atoms. In (a) 
and (b) the black dashed lines indicate the unit cell of each 
surface, (a-2) and (b-2) represent the triangular lattice used 
for the tight-binding models. In (a-2) the blue balls denote 
the Tl atoms, while in (b-2) each purple ball denotes a set 
of three Bi atoms. The red dashed lines indicate the mirror 
planes {xz plane). 


with respect to the xz plane Mxz [see Fig. [T][a-2)], the 
Hamiltonian is represented as 

H = ^ ^ cj[t + iXxip' X dij^z , (1) 

where a is the vector of the Pauli matrices, dij is the 
vector from site i to j, and Ci is an annihilation operator 
for an electron at the i-th site in the triangular lattice. 
^ij takes the values ±1 depending on the hopping direc¬ 
tions: ^ij is ±1 for 0 = (4nT l)7r/6 (n: integer), where 0 
is an angle between the hopping vector and the -\-x direc¬ 
tion. The second term of the Hamiltonian involves the in 
plane spin perpendicular to the hopping direction. The 
third term involves the spin perpendicular to the plane 
(Ijz), and such term does not exist in Rashba systems. It 
causes non-Rashba spin splitting with spin polarization 
perpendicular to the crystal surface at the A, i?' points, 
as we see in the following. 

The Hamiltonian m is rewritten into a matrix form 


H[k) — do diCTx + d20'y + d^cTz, (2) 

where 

do = t(2 cos2Y + 4cos A cos Y), (3) 

= Aa,(2 sin2Y-|-2 cos A sin Y), (4) 

^2 = — Aa,(2\/3sin A cos Y), (5) 

= A^(2 sin2Y — 4 cos A sin Y), (6) 

where A = ^kxa^ Y = \kya^ and a is the lattice con¬ 
stant. Its eigenvalue and spin direction are given by 

E{k) =do+ v^dj + dl+dl (7) 

(s) = rid, (8) 

where 7^ = ±1, d is a unit vector along d = (di,d 2 ,d 3 ). 


The surface Brillouin zone is a hexagon, with its corners 
A(0,1^), A'(0, —|j). As an example, we numerically 
calculate the band structure for t = 1, = 0.1. At A 

and K' points di = d 2 =0, and the Hamiltonian reduces 
to do + dscr^. Hence the energy bands around the A, A' 
points have non-Rashba splitting, with the spin direc¬ 
tions perpendicular to the crystal surface [Fig. [2][a)]. On 
the other hand, at M(^^=^, 0), M'(— 0) and r(0, 0), d 
becomes zero, and the two bands are degenerate. Around 
M, M' and F points, the vector d forms a vortex, and 
consequently the bands have Rashba splitting (Fig.[2]^b)). 
These features of our results qualitatively agree with the 
experimental results for the Tl/Si(lll) surface [T5|. The 
spin distribution for the present model in the first Bril¬ 
louin zone is shown in Fig. He). This is in agreement 
with the above considerations. 

Here, we note a previous work on theoretical calcula¬ 
tion for the Tl/Si crystal by a different method [2^. It 
is based on a four-band effective model on the honey¬ 
comb lattice representing the Tl atoms and Si atoms in 
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the first layer. The calculated spin configuration qualita¬ 
tively agrees with our results. Compared with the four- 
band model Hamiltonian, our two-band model is simpler 
and useful for analytic calculations and investigations of 
new phenomena. 


(a) (b) 



FIG. 2 . (Color online) Energy bands for effective surface 
Hamiltonian of Tl/Si. Band structures are plotted along (a) 
the K-T-K' and (b) the M-T-M' directions. Spins at K and 
K' are shown as blue and red arrows in (a), (c) Spin config¬ 
uration of the upper band (77 = 1) in the first Brillouin zone. 
Arrows represent the in-plane spin direction ((&), {Sy)) and 
the color represents the out-of-plane spin component {Sz)- 
The black hexagon represents the surface Brillouin zone. The 
parameters are set as t = 1 , Acc,;^ = 0 . 1 . 


B. /3-Bi/Si(lll) surface 


Next, we take the same procedure for the Bi/Si crystal 
surface. We construct an effective tight-binding Hamil¬ 
tonian on the triangular lattice, representing the states 
close to the Fermi energy. Symmetries of the Bi/Si crys¬ 
tal to be considered are the Csz symmetry, the time- 
reversal symmetry, and the mirror symmetry M^z [see 
Fig. [T][b-2)]. The symmetry properties might look the 
same with Tl/Si(lll) but it is not the case. As can be 
seen from Fig. [Ha-2) and (b-2), the relative positions 
of mirror planes in the triangular lattice are different, 
and it brings about different restrictions for the effective 
Hamiltonian. The resulting effective model with only the 


nearest-neighbor hoppings is written as 

H='^Cj[t + iXyia X (9) 

(*d) 

This Hamiltonian does not contain a term similar to the 
third term in Eq. ©• The Hamiltonian matrix can be 
written as 


H[k) = do dicjx + d2(7y^ (10) 

where 

do = t(2cos2X + 4cosX cosT), (11) 

di = A^(2a/3cosX sinF), (12) 

(i 2 = —Ay(2 sinX cosF + 2 sin2X), (13) 


where X = \kxa^ Y = ^kya^ and a is the lattice con¬ 
stant. As an example, we calculate the band structure for 
t = 1 and A^ = 0.1, shown in Fig. [3l Around the X, M, 
and F points, the vector (^ 1 ,^ 2 ) forms a vortex, and as a 
result the energy bands [Figs.[3][a) and (b)] have Rashba 
splitting around X, M, and f points. The spin distri¬ 
bution is in plane, and is shown in Fig. life). Usually, 
Rashba splitting appears around time-reversal-invariant 
momenta, where the spin degeneracy comes from the 
Kramers theorem. Nevertheless, in the present case, 
there is a Rashba splitting at the K point which is not in¬ 
variant under time reversal. It is called peculiar Rashba 
splitting [ 1 ^ which is different from the conventional 
Rashba splitting, because this peculiar Rashba splitting 
comes from the C^z symmetry at the K points which are 
not time-reversal invariant. These results qualitatively 
agree with the experimental results for /3-Bi/Si(lll) sur¬ 
face in Ref. M- 


III. CURRENT-INDUCED SPIN 
POLARIZATION IN THE TL/SI SURFACE 

As is already known, in systems with Rashba SOI, 
a charge current induces the spin polarization, be¬ 
cause an off-equilibrium electron distribution necessar¬ 
ily induces spin imbalance. It was theoretically pro¬ 
posed in Refs. [2ll-[^ and experimentally observed in 
Refs. [ 23 , [ill. For the Rashba systems the total spin 
polarization is in-plane, and is perpendicular to the cur¬ 
rent in the plane. Here we consider this current-induced 
spin polarization in systems with non-Rashba-type spin 
splitting, such as in the Tl/Si model. Because the third 
term in the Tl/Si model in the Hamiltonian ([1]) gives rise 
to the spins perpendicular to the crystal surface, we can 
naively expect the current-induced spin polarization to 
have an out-of-plane component. 

We numerically calculate the induced spin by the 
Boltzmann equation [ 2 ^. We assume that the model 
is uniform and that the impurity concentration is very 
small. In principle, the Boltzmann semiclassical equa¬ 
tion can be obtained from the Wigner quantum kinetic 



















4 



FIG. 3. (Color online) Energy bands and spins for effective 
surface Hamiltonian of the Bi/Si. Band structures are plot¬ 
ted along (a) the K-T-K' and (b) the M-f-M' directions, (c) 
Spin polarization for the upper band (77 = 1). Each arrow 
represents the in-plane spin direction ((&), {Sy)). The black 
hexagon represents the surface Brillouin zone. The parame¬ 
ters are set as t = 1 , A^, = 0 . 1 . 


equation as a semiclassical limit, and the classical distri¬ 
bution function is written as /cr,cr'(G which depends 
on spin indices. In this paper, our goal is to qualita¬ 
tively discuss various aspects of the non-Rashba SOC, 
and not to quantitatively calculate physical quantities in 
specific materials. Therefore, for simplicity of calcula¬ 
tion, we assume the classical distribution function to be 
spin independent, /(r, k, t). It means that the relaxation 
time is assumed to be independent of spin. In addition, 
the classical distribution function / does not depend on 
r, because we assume this system to be uniform. The 
Boltzmann equation becomes 


o = ^-v f+(^I 

dt ^\dt 


(14) 


collisions 


where the last term is the temporal change rate by col¬ 
lisions with impurities. The equation of motion for ^ 
in the absence of an external magnetic field is given 
by ^ = —eE. We approximate the collisions term by 
relaxation-time approximation: 


m - /o(fc) 

/ collisions '^t{k^ 


(15) 


where /o(A)) is the Fermi distribution function, and Tt{k) 


is the relaxation time. Here, we assumed that the relax¬ 
ation is due to impurity scattering, and is independent 
of spin, in order to illustrate the role of the non-Rashba 
SOC in the current-induced spin polarization in a simpler 
fashion. In particular, we neglect scattering by phonons 
which may lead to anomalous Cherenkov effect due to 
the interaction between electrons and the crystal lattice 
with SOC [13, because we assume the electric field not 
too strong so that the electrons are subsonic, and, thus, 
the scattering rate on impurities may dominate over the 
Cherenkov dissipation. By taking an approximation of 
replacing / by /o in the first term of Eq. ([H]) , the distri¬ 
bution function under an electric field is written as 

f{kMo{k) - ert{k)E • V^fo{k), (16) 

The spin polarization is written as 

{Si) = \j J^i^i) [-eTt{k)E ■ Vg/o(fc)) , (17) 

where {si) = r]di is the spin expectation value. 

The numerical results are plotted in Fig. 01 where we 
take t = 1 and = A;^ = 0.1 and set to be a con¬ 
stant r. Here we note that the nonzero spin polariza¬ 
tion in Fig. 0] arises from the spin expectation values for 
spin-split bands, because the relaxation time is set as 
a constant, and the spin dependence is considered only 
through the spin expectation value. It is also the case 
for the conventional spin-orbit-coupled systems such as 
Rashba systems |2ll-[^ . The electric field is taken along 
the X direction in Fig. 01^a) and along the y direction 
in Fig. 0|^b). The result is shown as the dimensionless 
spin polarization which is the spin polarization divided 
by = x,7/), as a function of the Fermi energy 

Ef. We see that for both cases the induced spins are 
along the direction which is the direction rotated by —90° 
from the electric field, which follows from symmetry as 
we show later. From these results, in spite of the presence 
of the third term in Eq. (|2]) having the out-of-plane spins, 
the current does not induce out-of-plane spin polarization 
within linear response. It means that the contributions 
from the K and K' points cancel each other. 

The results in Fig. 0] have three characteristic ener¬ 
gies: E^p = —3.5 where the spin polarization arises, 
E^P = —2.3 for the kink, and eJ"^ = —2 for the 
peak. From the energy band [Fig. [2], E^p = —3.5 and 
eP = — 2.3 correspond to band bottoms of the spin-split 
bands near the the K point (Fig. [2][a)). At Ep = — 2, 
there are band bottoms of spin-split bands at the M point 
(Fig. Efb)). To see the reasons for the characteristic val¬ 
ues of the Fermi energy Ef in Fig.[2l we calculate the con¬ 
tribution of the induced spin polarization, coming from 
the states near the K and that from the M points. The 
results for the spin polarization are shown in Fig. [5l We 
take the numerical calculation range of Fermi energy set 
as —2 > Ef^ because the approximate Hamiltonian is 
effective only in the vicinity of each point. 
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(a) 


(b) 



FIG. 4. (Color online) Spin polarization induced by the elec¬ 
tric field along the (a) x direction and (b) y direction. Spin 
polarization is shown in a unit of 


The contribution from the states near M point 
[Fig. [5][a-l)] has a peak at Ef = —2. Thus, the peak 
of the spin polarization in Fig. |l]is due to the contribu¬ 
tion from the states close to the M points. The strong 
peak due to the states close to the M point might be 
attributed to a relatively fiat dispersion around the M 
point, giving rise to a high density of states. On the 
other hand, the contribution from the K point (Fig.[5][a- 
2)), has a kink at Ef = —2.3. The kink at E^p is due to 
the contribution from the states close to the K points. 

(a-1) (a-2) 



FIG. 5. (Color online) Spin polarization induced by the cur¬ 
rent for the states around (a-1) M and (a-2) K in response 
to Ex- Spin polarization is shown in a unit of 

To see the reason for the absence of out-of-plane spin 
polarization in Fig. 01 we express the linear response of 
the spin polarization to the current as 

{Si) = (pijJj^ (18) 

where (pij represents response coefficients and Jj {j = 
x^y) is the current. By imposing crystal symmetries to 
Eq. (fT8|) . we obtain (pxy = —<Pyx and (pij = 0 otherwise. 
Thus, it is written as 

{S)=c^{Jxz), (19) 

where 0 is a real constant. From this calculation, we 
find that the spin polarization is always in-plane, and is 
always perpendicular to the in-plane current. It agrees 
with our numerical calculations and with the calculation 
by Liu et al. [2^ with four-band model. Therefore, in the 
present case the crystal symmetries prohibit the out-of¬ 
plane spin polarization, despite the nonzero out-of-plane 
spins for the states near K and K' points. 


IV. PERSISTENT SPIN HELIX IN THE TL/SI 


In the model of Tl/Si, the spin splitting around the 
K and K' points brings about an out-of-plane spin tex¬ 
ture around these points, similar to thin films of transi¬ 
tion metal dichalcogenides In such cases, nest¬ 

ing between the two pockets at K and K' brings about 
an anomalously long lifetime of a spin excitation at the 
nesting wave vector, as proposed in the two-dimensional 
electron gas with equal size of the Rashba and Dressel- 
haus spin-orbit couplings [12|. The key ingredients of 
this novel phenomenon are the two nested Fermi surfaces 
by the magic shifting vector having the opposite spins. 
In addition, the enhanced spin lifetime has been con¬ 
firmed experimentally in semiconductor heterostructures 


We apply this theory to the Tl/Si model. We assume 
that the Fermi energy is controlled by doping, so that 
there are small electron pockets around K and K'. Then 
the system has the Fermi surfaces shifted by the magic 
shifting vector Q = [0, —87r/(3a)] from the K point to 
the K' point [Fig. [6][a)] and the spins have the opposite 
direction perpendicular to the xy plane at the K and K' 
points [Fig. Efa)]. For example, let us take a state with 
the spin along the x-direction at x = 0, described by 



where the ko (|fco| ^ tt/o) is a wave vector. It is a su¬ 
perposition of two eigenstates at K and at K'^ which are 
degenerate. The spin expectation value at arbitrary x is 
written as 

{Sx) = = cos (q ■ , ( 21 ) 

(5'y) = f(V’k»/l^) =-sin(Q-f) , (22) 

(5'^) = f(V’k^lV’) = 0. (23) 

Hence, the rotation angle of spin moving from x = 0 to 

X p 0 is Q 'X. Therefore, the rotation angles of spin from 
X = 0 to X = di and 0*2 are 

4:71 ->■ Att 

Q-a-{ = —, Q-a 2 = - — . (24) 

From this calculation we obtain Fig. [6l[b) as a spin tex¬ 
ture on the T1 atoms forming a triangular lattice. It 
is a spin helix with the wave vector Q. As is similar 
to the spin helix in Ref. [T2|, it is expected to have en¬ 
hanced spin lifetime, because of the degeneracy of the 
states around K and those around K'. Therefore the spin 
helix in the Tl/Si model has in-plane spin texture with 
neighboring spins being different by 120°. This is dif¬ 
ferent from the spin helix in the 2D Rashba-Dresselhaus 
system , where the spin helix is perpendicular to the 
surface and a value of the magic shifting vector is much 
smaller than ours. To realize the persistent spin helix 
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in the Tl/Si system, hole doping [3l| leads to emergence 
of two nested hole pockets with opposite spins at the K 
and K' points. Once the spin texture is created, it will 
survive for a relatively long period. 



FIG. 6. (Color online) (a) The surface Brillouin zone (SBZ) 
in Tl/Si and the Fermi surfaces (blue circles) around K and 
K' points near Ef = —2. The red arrow indicates the magic 
shifting vector Q. (b) Schematic illustration of the spin helix. 
The orange arrows indicate the spin on the T1 atoms (blue 
balls). The directions of the red arrows represent the in-plane 
spin directions forming the spin helix. The vectors of al, oj 
(green arrows) indicate primitive translation vectors. 


V. BOUND STATE AT A JUNCTION 
BETWEEN TWO SURFACE REGIONS 

Next, we discuss bound states at a junction between 
two regions which have different signs of the SOI parame¬ 
ter, as in Fig.[7|^a). We consider the two models discussed 
in this paper, i.e., the Tl/Si system described by Eq. ([T]) 
without second term and the Bi/Si system described by 
Eq. (|9|). Eor the direction of the junction between two 
regions, we consider two cases shown in Eigs. El^b) and 
7(c). In the junction models, the SOI parameter (A;^ in 
Tl/Si and \y in Bi/Si) is set as +A and —A in the re¬ 
gions I and II, respectively. To extract interface states, 
we compare the results with 2D bulk models, where the 
SOI parameter is set as A for the whole system. 


A. Bi/Si junction model 

We calculated the band structures for two types of 
Bi/Si junction models shown in Eigs.[71^b) and 7(c). Eor 
Eig.[7Kc), we could not find bound states at the junction. 
In contrast, for Eig.[7|^b), we obtain the result in Eig.[8l^a), 
which shows that there are bound states (red lines) below 
all the 2D bulk states (gray area) in Eig.[8fa). The decay 
length depends on the SOI parameter; for example it is 
10.0 times the lattice spacing for the parameters jA^j = t. 
Additionally, we find that the energy difference between 
the bound states and the bulk band edge is larger for a 


(a) - 



Region I 
A>0 


Region II 
A<0 


FIG. 7. (Golor online) (a) Schematic illustration of the junc¬ 
tion model. The blue (I) [green (II)] area indicates the region 
of the Bi or T1 layer where the SOI parameter has a negative 
(positive) value. The yellow area indicates the Si substrate. 
In (b) and (c), two choices of junctions are shown as a top 
view. The primitive translation vectors (red arrow) are (b) 
(\/3a,0,0) and (c) (a, 0,0), respectively. The purple dashed 
line indicates the junction. 


larger value of the SOI parameter, as seen from the re¬ 
sults with the SOI parameter |A^| = 0.2t, |A^| = 0.5t, 
and |Ay| =t [Figs. [H][b)-8(d)). 

Furthermore, in Fig. [8][a), there are four bound states 
forming two Kramers pairs. The Kramers degeneracy ap¬ 
pears because the system restores spatial inversion sym¬ 
metry when the SOI parameter are set as —Xy{y < 0) = 
Xy{y >0). In contrast, when the SOI parameter is set 
as —Xy{y < 0) 7 ^ Ay(^ > 0), the degeneracy of the bound 
states are lifted at k points other than the M points. 

Additionally, spin directions of these bound states are 
almost along the z axis, with a small value of y com¬ 
ponent. This is expected from the crystal symmetries. 
Because this Bi/Si junction model [Fig. [7l[b)] has Myz 
and the time-reversal symmetries, the x component of 
spin necessarily vanishes. The out-of-plane spin direc¬ 
tions of the bound states are in sharp contrast to the 
in-plane spin directions of the bulk states. We note that 
these interesting bound states cannot be reproduced by 
continuum models as we discuss later. 

Next, we try to analytically calculate bound states for 
junction systems in order to compare with the above nu¬ 
merical results. First, we note that in the Bi/Si junction 
model, the bound state of Fig. [SJa) appears around the 
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FIG. 8. (Color online) Energy bands of the 2D bulk states 
and the junction model for Bi/Si with different SOI parameter 
such as (b) \Xy \ = 0 .2t, (c) \Xy \ = 0.5t, and (a), (d) |A^| = t. 
(a) The red lines indicate the states when the SOI parameter 
is set as in the blue region (I), 0 on the junction bound¬ 
ary and —t in the green region (II) (Fig. [7j). The gray areas 
indicate the bulk bands with the SOI parameter set as +t. 


M or M' points. These points are shown as green and 
blue arrows in Fig.[Hl^a), corresponding to the dotted lines 
in Fig.lTfa) and at these points the spin splitting is of the 
non-Rashba type. Therefore, we take the Hamiltonian 
from an expansion of Eq. around the M point: 

^2 

H (A:) — T Xy (^(Jxky T (^ykx) • 

We then consider a junction between two surface regions, 
described by the Hamiltonian (j25j) with the SOI param¬ 
eter changing along the y direction Xy = A(^). We tried 
the calculation of the bound states of this model, with 
two different choices of the function A(^): the step func¬ 
tion and the hyperbolic tangent function, as explained in 
the Appendix. As a result, the continuum model in the 
Bi/Si junction model cannot reproduce the bound states. 
The reason for this is yet to be clarified, and left as a fu¬ 
ture work. We note that for an interface between two 
regions of the one-dimensional nanowire with different 
directions of the SOI vector, bound states are found an¬ 
alytically [HI , whereas the system considered is different 
from ours. 

To realize this junction system, we note that the SOI 
parameter in this system originates from the broken 
inversion symmetry due to the surface, similar to the 
Rashba SOI. Therefore, realization of our junction model 
with different signs of SOI parameters would be similar 
to the case for the Rashba SOI. Here, we note an example 
of a related system, a non-centro-symmetric semiconduc¬ 
tor BiTel. In BiTel, the stacking order of atomic layers. 


Bi-Te-I or I-Te-Bi determines the sign of the Rashba pa¬ 
rameter in this system, and by changing the stacking 
order it is possible to achieve a different sign of Rashba 
SOI [HI- Although the present system is different from 
BiTel, it might be possible to achieve a junction between 
two regions with different signs of SOI parameters by 
controlling the atomic arrangement in each region. 


B. Tl/Si junction model 


Next, we numerically analyze two types of Tl/Si junc¬ 
tion models shown in Figs. ETb) and 7(c). There are 
no bound states in the Tl/Si junction model of Fig. [Il^c). 
Meanwhile, in the Tl/Si junction model of Fig.[71^b) there 
appear bound states (red lines) above all the 2D bulk 
states (gray area) in Fig. [9fa). The results for various 
values of the SOI parameter are shown in Figs.[9l^b)-9(d). 
From these results, the energy difference between the 
bound states and the bulk band edge is non-monotonic 
and it is maximized at about |A; 2 | = 1.73t [Fig. in 
our case. The decay length of the bound states also de¬ 
pends on the SOI parameter, and it is 38.5 times the 
lattice spacing, for the parameters |A; 2 | =t. 

As is similar to the Bi/Si junction in the previous sec¬ 
tion, the four bound states in Fig.[9l[a) form two Kramers 
pairs, stemming from the restored spatial inversion sym¬ 
metry when —Xz{y < 0) = Xz{y > 0) holds. This degen¬ 
eracy at k points other than the F point is lifted when 
—Xz{y < 0) ^ Xz{y > 0). We also note that the spin of 
the bound states is along the ^ axis because the effective 
Hamiltonian (Eq. ([1])) without the second term only has 
the az term. 

Next, we construct a continuum model for the present 
system and calculate bound states for the junction sys¬ 
tem in order to compare with the numerical results of 
the tight-binding models. The bound states for Fig.[7|[b) 
a,t kx = 0 [orange arrow at Fig. [9l[a)] should correspond 
to K, M, or F points [orange dotted line at Fig. [7|[b-2)], 
whereas in Fig. [7|[c) at kx = 0 corresponding to M and 
F points there are no bound states [orange dotted line 
at Fig. 0(0-2)]. Therefore, the bound states of Fig. [9Ka) 
are expected to come from the K point. For this reason, 
by expanding Eq. m around the K point we take the 
following continuum Hamiltonian: 

^ ^ 2m ^ ^ 

We take the same procedure for calculation of the bound 
states for this continuum model [Eq. (|26|) ] in the Ap¬ 
pendix. As a result, the continuum model cannot re¬ 
produce the bound states. Thus, we conclude that the 
bound states found in the junction system are unique to 
the lattice model. 
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in the Tl/Si model, the spins are shifted by 120° between 
the neighboring sites and are in-plane, which is different 
from the spin helix in semiconductor heterostructure. 

Finally, we numerically find bound states at the junc¬ 
tions between the two surface regions which have different 
signs of the SOI parameters in the Bi/Si system and in 
the Tl/Si system. For the junction shown in Fig. [Tfb) 
we found bound states at the junction. From the lattice 
model we find that the spins of the bound states in the 
Bi/Si model is out of plane, which is in contrast with 
in-plane spin distribution in the bulk. 

To summarize the whole results, surface states allow 
various types of non-Rashba-type SOI due to low sym¬ 
metries, and they are expected to lead to new spin prop¬ 
erties which are absent in conventional Rashba systems, 
and they have much room for future research. 
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VI. CONCLUSION 

Appendix 


In conclusion, we derive effective nearest-neighbor 
tight-binding Hamiltonians for surfaces of Tl/Si and 
Bi/Si by taking into account the crystal symmetries. The 
energy band of the Tl/Si model has non-Rashba splitting 
and the spin direction for each band is perpendicular to 
the crystal surface at the K point. The energy band of 
the Bi/Si model has peculiar Rashba splitting at the K 
point, because the Rashba splitting occurs around the K 
points which are not time-reversal invariant. Our results 
of the model calculation qualitatively agree with the ex¬ 
perimental results. In addition, differences between the 
Tl/Si model and the Bi/Si model lie in the out-of-plane 
spin at the K points in the Tl/Si model. This differ¬ 
ence originates from the non-Rashba term involving the 
out-of-plane spins in Eq. (HD. 

Additionally, as one of our theoretical explorations to¬ 
wards novel spin properties, we calculate current-induced 
spin polarization for the Tl/Si model. Although the third 
term in the Tl/Si model in the Hamiltonian (Eq. ([T))) has 
a spin component perpendicular to the crystal surface, 
the out-of-plane component of current-induced spin po¬ 
larization within linear response is zero, as prohibited by 
the crystal symmetries. As a result the current-induced 
spin polarization of the non-Rashba model is in-plane, 
an d perpendicular to the current, which is qualitatively 
similar to that of Rashba model. 

Furthermore, we apply the theory of persistent spin 
helix to the Tl/Si model. We find that in the spin helix 


In this appendix, we analytically calculate the bound 
states for the Bi/Si junction model of Fig.[71^b). First, we 
consider a junction between two surface regions described 
by the Hamiltonian (|25]) with different signs of the SOI 
parameter, \y{y) ±Ao for y Too. Several choices of 
\y{y) are possible. First, we take \y{y) to have the form 
of the step function along the y direction \y{y) = XoO{y), 
where 0{y) is step function. The system is set to be 
infinite along the x direction. We should replace ky by 
—idy because there is no translational symmetry along 
the y axis. Furthermore, we replace Xy{y)[(Tx{—idy)] by 
—i^cFx[Xy{y)dy + dy\y{y)] to preserve Hermiticity of the 
Hamiltonian. 

Along the y direction, we call the two regions I (^ > 0) 
and n {y < 0). From the Hamiltonian rewritten from 
Eq. (|25]) . boundary conditions are written as 


V^7(0)=V^/7(0), (A.I) 

-— {9y'0/(O) - 9y'0/7(O)} = -icra,Ao'07(O). (A.2) 


Meanwhile, when the wave function of the bound states 
decays exponentially as the eigen energy is 


^± = ^T-4)±Ao 


G 

2m 


± AoC±, 



(A.3) 
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where C± = y k‘^ — Here, because the values of 

Xy{y) for the regions I and 11 have the same size with 
different signs, the decay of the wave functions in region 
n is characterized by /^±, as is the same with region 1. 
Therefore, Ei± = En^ follows, and the bound states of 
two regions are written as 

xPi=a ^ ^-K+v+ik^x 

+ /3 ^ (A.4) 

xl)u =a' ^ ^>^+v+ik:,x 

+ /3' <«;- + fc:r) ^ g«_y+ifc.a.^ (^.5) 

where a, /3, a' and jS' are constants. We then impose 
the boundary conditions to derive a, /3, a' and /3' of Eqs. 
()A.4p and (|A.5p . As a result, we found that there are no 
such values which satisfy the boundary conditions. 

Second, we try with the SOC parameter Xy{y) to have 
the form Xy{y) = Aq tanh(^), where a represents a width 
of the domain wall. The results of the Bi/Si junction 
model are drawn as red lines and the bulk band is drawn 
as a gray area in Fig. flOfa). According to Fig. fTOfa). 
there are no states below all 2D bulk states around M 
point (kx = 0). Altogether, we cannot reproduce the 
bound states by this continuum model. 


Next, we analytically calculate the bound states for 
Tl/Si junction model of Fig.[Hl^b). We consider a junction 
between two surface regions described by Hamiltonian 
Eq. (|26|) with different signs of the SOI parameter along 
the y direction. We try with the SOI parameter A;^ to 
have the form A;^ = Aotanh(|). We take the calculation 
procedure used for the Tl/Si junction model and obtain 
the energy band in Fig. [TOlbb According to Fig. [TOlbb 
there are no states above all 2D bulk states around K 
point {kx = 0). Altogether, we cannot reproduce the 
bound states by this continuum model. 



(b 
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FIG. 10. (Color online) Energy band of the 2 D bulk states 
and the junction model for (a) Bi/Si and (b) Tl/Si. The 
red lines indicate the states when the SOI parameter is set 
as Aotanh(^). The gray area indicates 2D bulk states when 
the SOI parameter is set as Ao in all the regions. The other 
parameters are set as Ao = 40, m = 0.5, and a = 0 . 001 . 
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